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Abstract 

We analyze the stochastic response of a finite set of globally coupled noisy bistable units driven 
by rather weak time-periodic forces. We focus on the stochastic resonance and phase frequency 
synchronization of the collective variable, defined as the arithmetic mean of the variable character- 
izing each element of the array. For single unit systems, stochastic resonance can be understood 
with the powerful tools of linear response theory. Proper noise induced phase frequency synchro- 
nization for a single unit system in this linear response regime does not exist. For coupled arrays, 
our numerical simulations indicate an enhancement of the stochastic resonance effects leading to 
gains larger than unity as well as genuine phase frequency synchronization. The non-monotonicity 
of the response with the strength of the coupling strength is investigated. Comparison with simpli- 
fying schemes proposed in the literature to describe the random response of the collective variable 
is carried out. 
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I. INTRODUCTION 



The response of noisy nonlinear systems driven by weak external driving terms is fre- 
quently analyzed with the powerful tools of linear response theory (LRT) l|, |2|. The am- 
plification of the amplitude of the average response with the noise strength as well as the 
non-monotonic behaviors of the amplitude and of the output signal-to-noise ratio are man- 
ifestations of stochastic resonance (SR) [3]. These behaviors have been indeed rationalized 
using the ideas of LRT [4]. Within the linear regime described by LRT, the response of the 
non-linear system is necessarily very noisy. The output noise level is much higher and the 
decay time of the fluctuations much longer than those present in a linear system subject to 
the same noise and driving terms. Consequently, the signal-to-noise ratio of the nonlinear 
system in that regime is very small, even at its peak value, and the gain does not exceed 
unity [5|]. If one is interested in obtaining large average amplitudes and large signal-to- noise 
ratio, one must find a way to control the output fluctuations and this is only possible in a 
nonlinear regime. 

A possible way to reduce the output fluctuation levels in noise induced SR is to use 
external driving terms with amplitudes so large that, even though they are still subthreshold, 
they alter significantly the potential relief of the dynamics and render invalid the LRT 



assumptions [6]. Anotherpossibi 



nonlinear oscillators. In 
systems were studied. 
More recently 
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ity is to concentrate on the global response of a set of 



10| the enhancement of SR effects in arrays of bistable 
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12l |. we have analyzed the collective response of a finite set of globally 
coupled bistable systems. We demonstrated that SR is indeed much enhanced with respect 
to SR in single bistable units. Indeed, gains larger than unity were observed for subthreshold 
sinusoidal driving forces. Those findings indicated that the arrays were indeed operating in 
nonlinear regimes. 

Another aspect of the stochastic response of the system refers to its noise induced syn- 
chronization with the driving term. By a suitable definition of a phase associated to the 
random output process, one can define an output average phase frequency and a phase dis- 
persion. The matching of the phase frequency with the driving frequency for a range of noise 
values is what is termed noise induced frequency synchronization. It has been analyzed with 
a variety of analytical and numerical procedures Q, [isl, [l^ , even in circumstances where 



2 



quantum tunneling effects are relevant [17|]. In [18|] noise induced phase synchronization of 
the collective variable for wide ranges of noise values were also observed when the array was 
driven by subthreshold inputs. 

The results reported in the above mentioned work were obtained with driving terms that, 
even though they were unable to produce the cited effects in the absence of noise, they were 
large enough to invalidate LRT when applied to a single bistable system. The question we 
address in the present work is whether the enhancement of SR and synchronization effects in 
the collective response of finite arrays of bistable units still persist when the driving external 
force is rather weak. By weak, we will mean here that: i) the SR effects induced by the 
external driving in a single bistable unit can be satisfactorily described, even at a quantitative 
level, by the LRT approximation; and ii) noise induced phase frequency synchronization in 
a single bistable unit does not properly exist. 

In the next Section, we introduce the model system and define the relevant quantities 
that characterize the phenomenon of SR and phase synchronization. In Section UTTl we show 
the results obtained by numerically solving the dynamical stochastic equations. Compar- 
ison with the predictions of recently formulated approximate descriptions of the collective 



dynamics 
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201 1 is carried out in Section IV. The last section concludes with some remarks. 



II. MODEL AND DEFINITIONS 

We consider a set of N identical subsystems, each of them characterized by a variable 
Xi(t) (1=1,. ■ ■ ,N) satisfying a stochastic evolution equation (in dimensionless form) of the 



type 
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6 ^ 

Xi = Xi-x^,+- J2{x, - Xi) + V2DUt) + F{t), (1) 

The external driving force is periodic in t, F{t) = F{t + T). The term C,i{t) represents a 
white noise with zero average and {C,iit)C,j{s)) = 6ij6(t — s). We define a collective variable, 
S{t), 

^W = 4E^.W' (2) 
j 

and concentrate on its long time response when the system size, A^, is kept finite and the 
amplitude of the driving term is weak (in the sense that when the driving force acts on a 
single isolated unit, SR is well described by linear response theory). We define the one-time 
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correlation function 



L{t) = ^ r dt{S{t)Sit + T)),, (3) 



TJo 

The notation (...) indicates an average over the noise reahzations and the subindex * indi- 
cates the long time limit of the noise average, i. e., its value after waiting for t large enough 



that transients have died out. As shown in our previous work, 22| we have that 



L(r) = Lcoh(r) + Lincoh(r) (4) 

where the coherent part, Lco/i(t), is periodic in r with the period of the driving force, while 
the incoherent part, Lincoh{T) arising from the fluctuations of the output S{t) around its 
average value, decays to zero as r increases. The output signal-to-noise ratio. Rout is 

Rout = hm — ~ — — = ^ — , (5) 



+0+ 



^incoh 



(n) Lincoh(^) 



where Q is the fundamental frequency of the driving force F{t), Lcoh(^) is the corresponding 
Fourier coefficient in the Fourier series expansion of Lcoh(T), and Lincoh(^) is the Fourier 
transform at frequency Q of Li^icohi'r)- 

For a set of coupled linear oscillators driven by an external driving force F{t) and 
subject to the noise terms ^i(t) as in Eq. ([T]), the SNR of the corresponding collective 
process, R^^t^ coincides with that of the random process formed by the arithmetic mean of 
the individual noise terms ^j(t) plus the deterministic driving force F(t), namely, F(t) + ^(t) 
with ^(t) = A^~^ Y^iLi ^j(^)- The process ^{t) is a Gaussian white noise of effective strength 
D/N. Then, it is easy to prove that 

(^) _ 2A'N[1 - cos(7rr)] 



Thus, for our nonlinear case, it seems convenient to analyze the SR gain, G, defined as 



2^ 



R, 



G^^, (7) 

which compares the SNR of a non-linear system with that of a linear system subject to the 
same stochastic and deterministic forces. 

In the case of noninteracting units = 0), the SNR of the collective output is times 



larger than that of a isolated unit driven by the same force. Nonetheless, as discussed in 
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the gain associated to the collective output is just the same as the one of a single, isolated. 



unit. Thus, for the weak forces that we are considering here, we expect that the collective 
gain will not exceed unity. As seen below, our numerical results will confirm that result. 

Another aspect of the response is the noise induced phase frequency synchronization. We 
note that for low noise strengths and driving forces with sufficiently large periods, a random 
trajectory of S{t) contains essentially small fluctuations around two values (attractors) and 
random, sporadic transitions between them. For each realization of the noise term, we then 
introduce a random phase process, (f){t), associated to the stochastic variable S{t) as follows. 
We refer to a "jump" of S(t) along a trajectory, when a very large fluctuation takes the 
S{t) trajectory from a value near an attractor to a value in the neighborhood of the other 
attractor. We count N^^'^t), the number of jumps in the a trajectory of the process S{t) 
within the interval (0, t] . A trajectory of the phase process is then constructed as 



=7rArM(t), 

so that (f){t) increases by 2n after every two consecutive jumps. 
The first two moments of the phase process are estimated as 

1 M 
a=l 



(8) 



(9) 



v{t) 



M 

M 



- m? 



«=i 



M 
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Lo=l 



(10) 



where, Ai is the number of generated random trajectories. 

The instantaneous phase frequency is easily determined from the time derivative of {4>{t)). 
After a sufficiently large number of periods of the driving force, n, the system forgets its ini- 
tial preparation, but the instantaneous phase frequency is still a function of time. Then, we 



define a cycle average phase frequency fl 



over a period of the external driving. 
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by averaging the instantaneous phase frequency 
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_ ^ ^ .(n+i)T dim ^ (0[(n + l)r])-(0(nT)) 

T JnT dt T 



(11) 



Similarly, the cycle average phase diffusion coefficient is evaluated from the instantaneous 



slope of the variance v{t) as [l7|, [23 1 



1 

T JnT at 



V [{n + 1)T] - V {nT) 
T 



(12) 



In previous works, approximate analytical expressions for these two quantities have been 



derived for the = 1 problem in the classical 
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26 1 and quantum cases [jj 



Those expressions can not be applied to the collective variable of an iV-particle problem, 

III. NUMERICAL RESULTS 

In general, nonlinearities preclude exact analytical solutions of Eqs. ([T]). We will use 
numerical simulations to obtain useful information about the stochastic process S{t). In the 
asymptotic limit ^ oo, Desai and Zwanzig showed that the statistical properties of the 
model could be analyzed in terms of a nonlinear Fokker-Planck equation which allows the 
coexistence of several stable probability distributions for some ranges of noise strengths and 
9 < 1. In the same asymptotic limit, we analyzed a few years ago the stochastic resonant 
behavior of the first moment, (S'(t))*, when the system is driven by a time dependent sinu- 
soidal force, using a combination of analytical and numerical procedures [9i]. In particular, 
for very weak input amplitudes, a linear response theory analysis showed that a huge am- 
plification in the amplitude of the average output {S{t))^ with respect to that of the driving 
force could be achieved. 

In this work we concentrate on situations where: i) the number of subunits is finite; ii) 
the amplitude of the driving force is rather weak. We have considered an external periodic 
rectangular driving, 

F{t) = (13) 

where n(t) = \2t/T\, \_z\ is the floor function of z, i.e., the greatest integer less than or 
equal to z. In other words, F{t) = A {F{t) = —A) if t G [nT/2, {n + l)T/2) with n even 
(odd). As detailed in Gj the numerical algorithm used to integrate the Langevin equations 



follows one of the schemes put forward by Greenside and Helfand [27 1. 

Let us first consider the case of independent subunits, 9 = 0. In Fig. [1] we display the 
behavior with respect to the noise strength, D, of the signal-to- noise ratio and the gain 
of the output variable S{t) for an array of 10 noninteracting bistable units, driven by a 
weak amplitude rectangular force {A = 0.1) with fundamental frequency Q = 0.01. For 
comparison purposes, in Fig. [2] we present the results for the signal-to-noise ratio and the 
gain of the response of a single bistable units operating under the action of the same driving 
term as in the previous figure Clearly, SR is manifested in the non-monotonic behavior of 
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FIG. 1: Signal-to-noise ratio, Rs (circles) and gain Gs (squares) of the collective variable for an 
array of iV = 10 non-coupled, {9 = 0), bistable units driven by rectangular inputs with Q = 0.01 
and A = 0.1. 

the signal-to- noise ratio with respect to the noise strength in both figures. But, as expected, 
the R values for the = 1 case are very small, 1/10 times the R values for the = 10 
case. On the other hand, the gain in both cases has the same values, not exceeding unity 
as required by the linear response theory. The expected amplification of the R values are in 
agreement with the predictions of the central limit theorem which indicates a reduction by 
an 1/A^ factor of the output fluctuations of the array with respect to those of a system with 
a single unit. 

The interactions between the units bring up changes in the collective response as depicted 
in Fig. [3l Here we have coupled the = 10 bistable elements with a weak coupling strength 
6 = 0.2. The existence of coupling increases substantially the non-monotonic behavior of R 
vs. D of the collective output relative to the uncoupled units case in Fig. [H Also, the peak 
value is reached at higher values of D. And, more importantly, the gain is clearly above 
unity for a wide range of noise values. This fact indicates that the array is operating in a 
nonlinear regime even though the driving amplitude A is rather small. 

The dependence on the interaction strength 6 of R and G is explicitly demonstrated in 
the next figure Fig. HI Here, we depict the peak values of the signal-to- noise ratio, -Rmax, 
and the gain Gmax on the coupling strength 6 for an array of A^ = 10 elements driven by 
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FIG. 2: Signal-to-noise ratio Rs (circles) and gain Gs (squares) for a single bistable unit driven 
by rectangular inputs with = 0.01 and A = 0.1. 




D 



FIG. 3: Signal-to-noise ratio Rs (circles) and gain Gs (squares) of the collective variable for an 
array of AT = 10 bistable units driven by rectangular inputs with $7 = 0.01 and A = 0.1 and 
coupling strength = 0.2. 

a rectangular input with parameters Q = 0.01 and A — 0.1. The non-monotonic behavior 
with 9 is clear. As the coupling strength increases from zero the SR effects become more 
pronounced until they reach a maximum at around 9 — 0.5. Increasing further the coupling 
strength leads to a decrease in the peak values. 
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FIG. 4: The peak values of the signal-to-noise ratio (Rmax) and gain (Gmax) of a = 10 array 
of bistable units vs. the coupling strength 9. The driving term is rectangular with = 0.01 and 
A = 0.1. 

The signal-to-noise ratio depends on the behavior of the correlation function. In the 
lower panel in Fig. 0, we depict the maximum value of the coherent part of the correlation 
function (Ccoh) on the coupling strength 6 for those noise values at which the amplification 
is maximal. In the upper panel, the dependence of the initial value of the corresponding 
incoherent part Cincoh(O) with 6 is shown. One observes that both quantities show non- 
monotonic behaviors with the coupling strength. Even though the peak in the lower panel 
and the minimum in the higher panel are obtained for 6 = 0.2, the largest signal-to-noise 
ratio is reached at the slightly larger coupling strength 6 = 0.5. This can be understood by 
noting that, as depicted in Fig. [5], the decay of the incoherent part of the correlation function 
at ^ = 0.5 is faster than at the other coupling values. Indeed, the key behavior is that of the 
incoherent part of the correlation function. Comparing Cincoh(^); for = and 6 = 0.2 we 
see that their initial values get smaller as 6 increases, while the decay time of the correlation 
function remains practically the same. On the other hand, as 6 increases to 0.5, the initial 
value also increases but the decay is much faster. Consequently, the denominator in the 
ratio of Eq. ([5]) decreases and the ratio reaches a maximum value. As the coupling constant 
is increased further, the initial values increase and therefore the values of R decrease. These 
features indicate that large enhancements of the SR quantifiers and large gain values are 
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FIG. 5: The maximum value of Ccoh{t) and the initial value of the incoherent part, Cincoh(O) for 
several values of 

achieved when the output fluctuations are small and fast decaying. These two combined 
facts can only be achieved when the system operates in a nonlinear regime. For a weak 
driving force this nonlinear regime is not possible with a single unit system or with an array 
of noninteracting units. 

Let us know turn our attention to the noise induced synchronization. In Fig. [8] the 
phase frequency (fiph) and the phase diffusion coefficient ( -Dph) obtained from numerical 
simulations are depicted. The driving force in all cases is a rectangular input with amplitude 
A = 0.1 and fundamental frequency Q = 0.01. It is clear that for this weak driving force 
and for a single unit system, the range of noise values at which the phase frequency matches 
the driving fundamental frequency is extremely narrow to consider that there is a proper 
synchronization. On the other hand, for = 10 coupled particles, the range of noise values 
leading to proper frequency matching and small phase diffusion coefficients is substantial. 
Optimum synchronization is obtained for 6 = 0.5 with a dip at the diffusion constant for 
D ^ 0.15. This noise value is slightly smaller than the one at which the signal-to-noise 
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FIG. 6: Coherent part of the correlation functions of the collective variable S{t) for = 10 and 
several values of 6. The driving term is rectangular with 0, = 0.01 and A = 0.1. The noise values 
used correspond to the values at which R reaches its peak for the different coupling strength. 
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FIG. 7: The phase frequency and the phase diffusion coefficient as functions of D. The driving 
term is rectangular with 0, = 0.01 and A = 0.1. 
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ratio reaches its maximum value for the same couphng constant {D ^ 0.19). Note that 
as the value of 6 increases above 0.5, the range of noise values for synchronization to take 
place reduces and it shifts to higher values of D. Perhaps, the most surprising result is that 
for uncoupled units, synchronization between the collective variable and the driving term is 
lost. This is indeed in sharp contrast with the results reported in [l^ for the same system 
with a driving force with a much larger amplitude A = 0.3. 

The range of noise values for stochastic resonance and noise induced phase synchroniza- 
tion do not have to coincide. Indeed, our results indicate that SR might be present at 
parameter values where noise induced phase synchronization does not exist. This is not 
surprising. Both effects are aspects of the stochastic response of the system to a driving 
agent. They are undoubtedly related, but they probe different aspects of that response. 



IV. COMPARISON WITH APPROXIMATE DESCRIPTIONS 



Due to the nonlinearity of the dynamical equations, the collective variable does not obey 
a single closed Langevin equation, but an infinite hierarchy of equations. Recently, for 
finite arrays, approximate descriptions of the dynamics based on truncations of the infinite 
hierarchy of fluctuating cumulants have been put forward by Pikovsky et al. [19] and by 
Cubero 2Q| . The Gaussian _approximation amounts to describe the dynamics of the collective 
variable by the equations |l9l |. 



2D 

S = S-S^- 3MS + J — e(t) + F{t) (14) 
= M- 3S^M - 3M^ -OM + D. (15) 

Even within the reduced Gaussian approximation, one has to rely on numerical treatments 
to obtain reliable information. The predictions of the Gaussian approximation (without 
invoking a slaving principle) are compared with those obtained from the full solution of the 
entire dynamics in the next figures. 

In Figs. [HI El M we show the results obtained for the noise dependence of the collective 
signal-to-noise ratio as given by the simulation of the full set of Langevin equations and by the 
Gaussian approximation. In Fig. [8] we consider the case of uncoupled arrays, while in Figs. 
[9] and [TO] we consider moderate (6* = 0.5) and strong {6 = 1.5) coupling cases respectively. 
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FIG. 8: Comparing the signal-to-noise ratio as obtained from the numerical simulation of the full 
set of equations and from the Gaussian approximation for = 10 uncoupled units( ^ = 0). Other 
parameter values: = 0.01 and A = 0.1. 

The Gaussian approximation seems to be more reliable as the value of the coupling strength 
increases, although it yields poor results around the noise values at which R shows its peak. 
Only for high values of the noise strength does the Gaussian approximation lead to results 
in good agreement with those obtained from the full set of equations. 



As discussed in [l9[, the Gaussian truncation combined with a slaving principle, allows 
for a simplified description of the dynamics of S{t) in terms of a Langevin equation in an 
effective double well potential and a noise term \f^i{t)- This proposal (or the effective 
potential proposed in [^l), has the attractive feature of reducing the dynamics for the col- 
lective variable to a single Langevin equation in a one- dimensional potential. Unfortunately, 
by contrast with the parameter values for 9 and D considered in those works, for the param- 
eter values used here the effective potential evaluated as indicated in the above mentioned 
references might not be bounded or it might not even exist. 

We have also compared the results obtained with the Gaussian truncation for the average 
phase frequency and the phase diffusion coefficient with respect to the ones obtained from 
the full set of Langevin equations. In Fig. [TT], we depict the results obtained for = 10 
uncoupled units {9 = 0) driven by a rectangular periodic force with A = 0.1 and fundamental 
frequency Q = 0.01. The Gaussian approximation predicts a perfect matching of the driving 
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FIG. 9: Comparing the signal-to- noise ratio as obtained from the numerical simulation of the full 
set of equations and from the Gaussian approximation for = 10 weakly {6 = 0.5) coupled units. 
Other parameter values: O = 0.01 and A = 0.1. 

frequency and the phase frequency for a wide range of noise strength values. These results 
are completely at variance with those obtained from the numerical solution of the full set of 
equations. Thus, the Gaussian approximation is not reliable for the 6' = case. 

On the other hand, as depicted in Fig. [121 the introduction of even a small coupling 
between the units drastically changes the picture. The Gaussian approximation results are 
in very good agreement with those obtained from the full set of equations. Indeed, as 
depicted in Figs. [T2|13|14|[T5l for coupled units the results of the Gaussian approximation 
are very good compared with those obtained from the full set of dynamical equations for 
some range of noise values. The range of noise values for phase synchronization depends on 
the coupling strength. As 9 is increased the range decreases. In all cases, synchronization 
starts disappearing as the noise strength becomes large. This is to be expected as large 
noise values might induce jumps over potential barriers quite independently of the driving 
force. Thus, the synchrony between noise induced jumps and the changes of sign of the 
external amplitude tends to disappear. The very same idea of the phase process introduced 
previously in Eq. ([8]) loses its meaning for large noise strengths 



14 
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FIG. 11: Comparing average phase frequency and phase diffusion as obtained from the numerical 
simulation of the full set of equations and from the Gaussian approximation for N = 10 and = 0). 
Other parameter values: O = 0.01 and A = 0.1. 



V. CONCLUDING REMARKS 



We have analyzed different aspects of the stochastic collective response of a finite array 
of globally coupled bistable units to a weak driving time-periodic force. We focus our anal- 
ysis on the phenomenon of stochastic resonance and noise induced phase synchronization. 
As demonstrated by our numerical results, both effects might be present in the collective 
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FIG. 13: Comparing average phase frequency and phase diffusion as obtained from tlie numerical 
simulation of the full set of equations and from the Gaussian approximation for = 10 and 
{9 = 0.5). Other parameter values: = 0.01 and A = 0.1. 

response as long as the units are not statistically independent, i. e., their coupling is not 
zero. 

There are two relevant facts: i) the gain of the collective variable might reach values 
greater than unity and, ii) the phase frequency might synchronize with the fundamental 
driving frequency in wide ranges of the noise strength. These two features clearly indicate 
that, for the weak driving forces considered here, the response of the system can not be 
analyzed with the tools of linear response theory. This is so, even though for a single unit 
subject to the same weak driver, LRT provides a very valuable tool to understand SR at 
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FIG. 15: Comparing average phase frequency and phase diffusion as obtained from the numerical 
simulation of the full set of equations and from the Gaussian approximation for = 10 and 
{9 = 1.5). Other parameter values: = 0.01 and A = 0.1. 

qualitative and even quantitative levels. Our calculations indicate that the failure of LRT 
is essentially due to the strong modification of the output fluctuations brought up by the 
external driving and the coupling between the elements of the array. The output fluctuations 
become much smaller and short lived than those present in a single unit system. 

We have also analyzed the non-monotonic behavior of the signal-to-noise ratio with the 
coupling strength when other parameter values (number of particles, amplitude and period 
of the driver) are kept constant. The largest signal-to-noise ratio occurs at an optimum 
value of the coupling strength. 
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Finally, we have also compared our numerical findings with those obtained by simplifying 
approximations that have been put forward in the literature. In particular, the closure 
of an infinite hierarchy of fiuctuating cumulant moments at the Gaussian level provides 
a rather good description of the simulation results at least for some range of parameter 
values. Unfortunately, further simplifications leading to effective one-dimensional Langevin 
dynamics for the collective variable, that seems to provide useful insight on the response for 
some regions of parameter space, become invalid in the parameter region considered in our 
work. 

Acknowledgment s 

Wc acknowledge the support of the Direccion General de Ensefianza Superior of Spain 
(BFM2005-02884) and the Junta de Andalucia. 



[1] R. Kubo, J. Phys. Soc. Jpn., 12, 570 (1957)) 

[2] Peter Hanggi and Harry Thomas, Phys. Rep. 88, 207 (1982)) 

[3] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, Rev. Mod. Phys., 70, 223 (1998). 
[4] M. I. Dykman, R. Mannella, P. V. E. McClintock, and N. G. Stocks, Phys. Rev. Lett. 65, 

2606 (1990); M. I. Dykman, R. Mannella, P. V. E. McClintock, and N. G. Stocks, Phys. Rev. 

Lett. 68, 2985 (1992) 

[5] Jesiis Casado-Pascual, Glaus Denk, Jose Gomez Ordonez, Manuel Morillo, and Peter Hanggi, 

Phys. Rev. E 67, 036109 (2003). 
[6] Jesiis Gasado-Pascual, Glaus Denk, Jose Gomez Ordonez, Manuel Morillo, and Peter Hanggi, 

Phys. Rev. E 68, 061104 (2003). 
[7] Peter Jung, Ulrich Behn, Eleni Pantazelou, and Prank Moss Phys. Rev. A, 46, R1709 (1992). 
[8] A. Bulsara and G. Schmera, Phys. Rev. E, 47, 3734 (1993). 

[9] M. Morillo, J. Gomez-Ordofiez, and J. M. Casado, Phys. Rev. E, 52, 316 (1995); J. M. Casado 
and M. Morillo, Phys. Rev. E, 52, 2088 (1995). 
[10] Alexander Neiman, Lutz Schimansky-Geier, and Prank Moss, Phys. Rev. E, 56, R9 (1997). 

[11] Jose M. Casado, Jose Gomez-Ordohez and M. Morillo, Phys. Rev. E 73, 011109 (2006). 



18 



[12] David Cubero, Jesus Casado-Pascual, Jose Gomez-Ordonez, Jose Manuel Casado and M. 
Morillo, Phys. Rev. E 75, 062102 (2007). 

[13] A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization: A Universal Concept in Non- 
linear Science (Cambridge University Press, Cambridge, 2001). 

[14] J. Kurths, S. Boccaletti, C. Grebogi, Y. C. Lai, Chaos 13, 126 (2003). 

[15] L. A. Vainstein and D. E. Vakman, Frequency Analysis in the Theory of Oscillations and 

Waves (Nauka, Moscow, 1983). 
[16] L. Schimansky-Geier, V. Anishchenko, and A. Neiman, in: N euro-informatics, edited by S. 

Gielen and F. Moss, Handbook of Biological Physics, Vol. 4 (Elsevier Science, new York, 2001); 

B. Lindner, J. Garcia-Ojalvo, A. Neiman, L. Schimansky-Geier, Phys. Rep. 392, 321 (2004). 
[17] Goychuk I. et al., Phys. Rev. Lett. 97, 210601 (2006). 210601. 
[18] J. M. Casado, J. Gomez- Ordohez and M. Morillo, EPL 79, 50002 (2007). 
[19] A. Pikovsky and A. Zaikin, M. A. de la Casa, Phys. Rev. Lett. 88, 050601 (2002) 
[20] David Cubero, Phys. Rev. E 77, 021112 (2008). 

[21] Rashmi C. Desai, and Robert Zwanzig, J. of Stat. Phys., 19, 1 (1978). 

[22] Jose M. Casado, Jose Gomez Ordonez, and Manuel Morillo, Phys. Rev. E 73, 011109 (2006). 
[23] J. Casado-Pascual, J. Gomez Ordonez, M. Morillo, J. Lehmann,I. Goychuk and P. Hanggi, 

Phys. Rev. E 71, 011101 2005. 
[24] J. A. Freund , A. B. Neiman and L. Schimansky-Geier , Europhys. Lett. 50, 8 (2000). 
[25] P. Talkner, Physica A 325, 124 (2003). 

[26] T. Prager and L. Schimansky-Geier , Phys. Rev. E 71, 031112 (2005). 
[27] H. S. Greenside and E. Helfand, Beh Syst. Tech. J., 60, 1927 (1981). 



19 



